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I. INTRODUCTION 



Non-equilibrium dynamics have become a very active field of research in the last few years 
in nearly all parts of physics. In condensed matter physics for example the description of 
the dynamics of non-equilibrium phase transitions plays an important role . Such phase 
transitions occur in ferromagnets, superfluids, and liquid crystals to name only a few. They 
are subjects of intensive studies, both theoretical and experimental. 

Also in cosmology some phenomena require the use of non-equilibrium technics. One 
example is the electroweak phase transition. The electroweak phase transition could lead 
to a possible model for baryogenesis. This problem is investigated by several groups, e.g., 
H [I]], using non-equilibrium methods. 

Another phenomenon in cosmology where non-equilibrium dynamics are important is 
the inflationary phase of the early universe which is studied intensively by different groups 
@> 0; 05 S B H3' H3' 03' EH' 0' Inflation refers to an epoch during the evolution of 



the universe in which it underwent an accelerated expansion phase. This would resolve 
some of the short comings of the Standard hot Big Bang model, e.g., the flatness problem, 
concerning the energy density of the universe and the horizon problem, related by the large 
scale smoothness of the universe, indicated by the Cosmic Microwave Background Radiation 
(CMBR). 

At lower energies heavy ion collisions are under consideration as non-equilibrium processes 



15| , |16| , |T7| , |18| , |19| , [20fl . In such heavy ion collisions a new state of matter could be reached 
if the short range impulsive forces between nucleons could be overcome and if squeezed 
nucleons would merge into each other. This new state should be a Quark Gluon Plasma 
(QGP), in which quarks and gluons, the fundamental constituents of matter, are no longer 
confined, but free to move around over a volume in which a high enough temperature and/or 
density reveals. Heavy ion collisions are studied experimentally at current and forthcoming 
accelerators, the Relativistic Heavy Ion Collider RHIC at Brookhaven and the Large Hadron 
Collider LHC at CERN. The occurring Quantum Chromodynamic (QCD) phase transition 
in these processes could be out of equilibrium and lead to formations of coherent condensates 
of low energy pions, so called Disoriented Chiral Condensates (DCC). Recent results reported 
from CERN-SPS |^TJ seem to indicate a strong evidence for the existence of a QGP in Pb-Pb 
collisions. 

The underlying theories of the above described phenomena are gauge theories. Therefore, 
it is of great importance to implement gauge field fluctuations if one wants to deal with more 
realistic models. Gauge theories out of equilibrium are still poorly explored in the literature. 
While the aspects of various gauges and of gauge invariance have been discussed extensively 
in perturbation theory and for equilibrium systems, the analysis of such systems out of 
equilibrium is still to be developed. The analysis presented here focus on different approaches 
for the description of non-equilibrium gauge theories and compare these approaches. 

In |22[ we have presented the theoretical framework and its numerical implementation for 
the SU(2) Higgs model in the 't Hooft-Feynman background gauge. Obviously, the question 
of gauge dependence and independence has to be explored. In [^] we investigated the gauge 
invariance of the one- loop effective action of a Higgs field in the SU(2)-Higgs model for a 
time independent problem, the bubble nucleation. An analogous analysis for a system out 
of equilibrium [E3 has shown that the fluctuation operator decomposes in a similar way to 



the one derived in p3|. We have found within this analysis a close relation between the 't 



Hooft-background gauge and the Coulomb gauge. 



2 



In this work we are now considering two other approaches in order to investigate gauge 
theories out of equilibrium and to get some insight how the behavior of the system is in- 
fluenced by the choice of gauge. Therefore, we study the Abelian Higgs model in a gauge 



invariant approach which was developed by Boyanovsky et al. |25j in order to study non- 
equilibrium aspects of the effective potential. We use their formalism in order to derive the 
equations of motion in the one loop approximation. We will indicate that the gauge invariant 
approach leads to IR-problems beyond the one loop approximation. As a second approach 
we examine the Coulomb gauge. Here we use a gauge condition in order to eliminate the 
unphysical degrees of freedom. We show how a proper renormalization can be implemented. 
Therefore, we use a scheme which was developed by us for the 4 theory and extended 



to fermionic and gauge systems p2| , |27| . Furthermore, we show how the gauge invariant 
approach and the Coulomb gauge are related to each other. 

In order to compare the behavior of the system for different gauges quantitatively, we 
reexamine the 't Hooft-Feynman gauge and implement the 't Hooft-Feynman gauge and 
the Coulomb gauge numerically. We examine the development of the zero mode and the 
different fluctuations. We also investigate the problem of energy conservation. 

The paper is organized as follows. In section 2 we start our discussion by outlining the 
main aspects of the gauge invariant approach of Boyanovsky et al. We derive the equation 
of motion for the zero mode and the different fluctuations and we determine the energy. In 
section 3 we present the equation of motion and the energy in the Coulomb gauge. The 
comparison of the two approaches is investigated in section 4. In section 5 we show how 
the renormalization can be done properly in order to prepare the equations for numerical 
implementations. Some numerical results for the Coulomb gauge and for comparison also in 
the 't Hooft-Feynman background gauge are presented in section 5. We conclude in section 
6 with a discussion of the influence of the different degrees of freedom and an outlook to 
more realistic applications of the method. 



II. A GAUGE INVARIANT APPROACH 
A. The Formalism 



We study the Abelian Higgs model in a gauge invariant formulation. The basic ideas for 
this description are developed in |1|, We give here a short overview about the derivation 



of the Hamiltonian. For details the reader is referred to the two papers mentioned above. 
The Lagrangian density for the Abelian Higgs model reads 

C = -\F^F, U + D^Dy - ±{<jt<j> _ , (2.1) 

D„ = 0„ + ieA^ . (2.2) 

We want to formulate a gauge invariant Hamiltonian. Therefore, we use as in fl], |25| the 
canonical formulation. We have to identify the canonical field variables and constrains. The 
canonical momenta conjugate to the scalar and vector fields are given by 

n° = , (2.3) 

IT = A i + VM° = -BP , (2.4) 

7r f = <p + ieA°0 , (2.5) 

7i = ft - ieA°ft . (2.6) 
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Therefore, the Hamiltonian reads 



H 



d 3 x \ ^fl • n f + vr+vr + (V0 - ieA<p) ■ (V0 t + ieA^) 



+~(Vx^ 



A 
4< 



v 2 ) 2 + A 



V • n - ie(vr0 - -n ] (p ] ) 



(2.7) 



We will quantize this system with Dirac's method [pqj . Therefore, we have to recognize the 
first class constraints (mutually vanishing Poisson brackets). Then the constraints become 
operators in the quantum theory and are imposed onto the physical states, thus defining 
the physical subspace of the Hilbert space and gauge invariant operators. We have two first 
class constrains 



and Gauss' law 



IT 

Q{x,t) -- 
P -- 



8C 

Ja® 



o, 



vv - p = 

ie (^(j)Tx — 0V 



(2.8) 

(2.9) 
(2.10) 



with p being the matter field charge density. We can now quantize the system by imposing 
the canonical equal time commutation relations 



[n°{x,t),A°(y,t)] = -i6(x-y), 
[W(x,t),A>(y,t)} = -i6«6(x-y) 
[^(x,t),(f)\y,t)} = -%8{x-y) , 
W(x,t),<f>(y,t)] = -%5{x-y) . 



(2.11) 
(2.12) 
(2.13) 
(2.14) 



In Dirac's formulation, physical operators are those that commute with the first class con- 
straints. Since U°(x,t) and Q(x,t) are generators of local gauge transformations, operators 
that commute with the first class constraints are gauge invariant [|T], ^3J. As shown in [JTJ] 
the fields and the conjugate momenta can be written in the following form 



$(£) 
Tl(x) 



(p{x) exp 
7r(x) exp 



ie d 3 yA(y) -V y G{y-x) 



d 3 yA(y)-V y G(y 



—ie 



x ) 



with G(y — x) the Coulomb Green's function that satisfies 

AG(y-x) = 5 3 (y-x) . 



(2.15) 
(2.16) 

(2.17) 



They are invariant under gauge transformations |T[|. The gauge field can be separated into 
transverse and longitudinal components 



A(x) = A L (x) + A T {x) 

V x A L {x) = , 

V ■ A T (x) = . 



(2.18) 
(2.19) 
(2.20) 
(2.21) 
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Since the fields Ax and $ and their canonical momenta commute with the constraints, they 
are gauge invariant. It is also possible to write the momentum canonical to the gauge field 
in longitudinal and transverse components 

n(f) = U L (x ) + U T (x ) , (2.22) 

where both components are gauge invariant. In jy, it is mentioned that in all matrix 
elements between gauge invariant states the longitudinal component can be replaced by the 
charge density 

u L (x) - le [$(m(y) - ^(y)n f (y)] = p ■ (2.23) 

Finally, the Hamiltonian becomes 

H = Jd 3 x |^I1 T ■ I1 T + n f n + (V$ - ieA T <$>) ■ (V$ f + ieA T &) +i(Vx A T ) 2 

+±(&<S>- v 2f+~Jd 3 x Jd 3 yp(x)G(x-y)p(y)} . (2.24) 

The features of this Hamiltonian are discussed at length in [[I], p5| . One of the striking 
points is the equivalence with the Hamiltonian in the Coulomb gauge. This similarity is 
not uncommon because, in the Coulomb gauge, only the physical degrees of freedom are 
taken into account. We focus our interest on the non-equilibrium aspects of the theory. In 
0, 25], the effective potential was derived and some aspects of non-equilibrium dynamics 
were discussed. But as there explain the effective potential is not suitable for non-equilibrium 
dynamics. We derive here the full non-equilibrium equations. We include not only the terms 
which are quadratic in the zero mode, but all terms up to second order which are relevant 
for the one loop approximation. We consider the derivative terms of the zero mode and 
of its conjugate momentum. This yields the wave function renormalization which is not 
considered in jl], |25| or has to be introduced by hand. We will also see that the formalism 
does not give a clear statement about the loop order which is included. When linearized, the 
equations become equivalent to those obtaining in the Coulomb gauge, which we consider in 
detail in section |I^. There, we show the correspondence of the Hamiltonian approach and 
the Coulomb gauge. As we will show, the inclusion of higher loop terms in the Hamiltonian 
approach leads to problems in the IR-region. 

First of all, we derive the equation of motion for the fields. Therefore, we separate the 
expectation value of <3? and of its canonical momentum into a mean value and a fluctuation 
part 

$(f,t) = <j>(t) + <p{x,t) , (2.25) 
U(x,t) = U(t) + ir{x,t) . (2.26) 

We also introduce real fields and canonical momenta as follows 

<P = + > ( 2 - 27 ) 

1 

Therefore, we find as a gauge invariant Hamiltonian 
1 



h = n 



2 n ' + U(<f>) 



Ti = —(n h -m v ) . (2.28) 



^2.29) 
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+ jd 3 x 



2 L 7"y,„\2 i ^ / . ,,2 ,,■ 7,2 , . 2 ■, > .2/. 2 



2 

+ |rf 3 x |d 3 y|-[^(f)-n^(f)]G'(f-y)[07r^(y)-n)^(y)] , 



with 



W) = t(0 2 -- 2 ) 2 . 



(2.30) 



B. Equation of Motion and Energy Density 



We are now able to derive the equations of motion for the zero mode and the fluctuation 
fields. In the Hamiltonian formalism we find the equation of motion by calculating the 
commutator between the Hamiltonian and the field. We obtain for the zero mode <p 



n 



i + 

2 



e 

and for the canonical momentum 



6 - Jd 3 x J d 3 y <p(x)G(x - y)<p(y 
d 3 x J d 3 y [ip(x)G(x - y)n v (y) + n v (x)G(x - y)<p 



(2.31) 



n 



-u 



(2.32) 



'(0) -f^Jd 3 x [e 2 d\ + X(h 2 + y? 2 ) + 2Xh 2 

2 

-^<pjd 3 x Jd 3 yir v (x)G(x-y)iT v (y) 

+ ^} U I d3x I d3y W*)<?(* - yMy) + <p{x)G{x - y)7r v (y)} . 

We also need the equations of motion for the three different quantum fluctuations. The 
first one is the transverse gauge field. We find an equation for the field a±i itself and for its 
canonical momentum. We can combine these two expressions to a second order differential 
equation for the gauge field: 



d±i = TT ±i , 7T±i = (A — e 

=>- d ±i = (A - e 2 4> 2 )a 1 _i . 



2 ±2 



(2.33) 



In the same way, we get the equation of motion for the real component of the Higgs fluctu- 
ation 



h 



A 



2 -v 2 ) 



h . 



(2.34) 



More difficulties arise for the Goldstone sector ip because the field and its canonical momen- 
tum are coupled via Green functions. We find for the field itself 

e 2 r 

(p = ir ip + — J d 3 x [(j)G xy ((f)Tx^ - U<f) + (07r v - Uip)G xy (f)\ , (2.35) 
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and for the momentum 



A<p - Xip{4> 2 



- Id 3 x 



U 2 G xyV - (11(f) + (pU)G xy 7i v 



[2.36) 



Now we transform the equations of motion for the zero mode and the fluctuations into 
Fourier space. Therefore, we expand the fluctuation fields in the following way 



i)(x, t) 



d 3 k 



— \a k U{k, f) e *-*-*"o* + a{U*(k, t ) e - its+lU]ot } 



[2.37) 



(2tt) 3 2lu 

with the usual commutator relations for the annihilation and creation operators ak,a\: 

[a h , a{,] = (2ir) 3 2uj 5 3 (k - k') . (2.38) 

The U(k, t) are the mode functions for the fluctuations depending on k. For convenience 
we omit the momentum dependence in the following. The Fourier transform for the Green 
function leads to a factor 1/k 2 . With these expansions, the equation of motion for the zero 
mode reads 



Ht) = n(t) 



l + e< 



d 3 k 



1 



(2tt) 3 2upok 2 
' rI '' ' %{t)m(t) + U^.(t)U\lt) 



(2tt) 3 2u p0 k 2 



[2.39) 



and for the canonical momentum 

ti(t) = -u'{4>) 



-4>(t) 

+U(t) 



d 3 k 
(2^)3 
d 3 k 



3A 



^2 



U ± (t)\ 2 + 7^-\U h (t)\ 2 + 



2u 



ha 



2uj 



l^)| 2 + 



2uj p0 k 2 



KM 



u^{t)u*Jt) + u;(t)u v {t) 



(2tt) 3 Auj pQ k 2 

where we have introduced the following frequencies 



(2.40) 



2 


= 'k 2 + X(<j ) l-v 2 ) 


(k 2 + e 




= P + X^l-v 2 ), 






= k 2 + m 2 w + e 2 ((j) 2 


-v 2 ), 




= k 2 + m 2 h + 3X(0 2 O 


~ v 2 ) , 


m 2 h 


o \ 2 2 2 

= 2Xv , m w = e 


2 

V . 



2 J,2\ 



2 2 
•Vo^aO 



(2.41) 
(2.42) 
(2.43) 
(2.44) 
(2.45) 

The index indicates the choice of t = 0. Notice, that we have included a factor two for the 
two transverse gauge freedoms in the equation of motion for IT. u p0 was introduced in [|TJ. 
By comparing these results with the zero mode equation in the gauge fixed theory which we 
have examined in [22] we find some analogies in the fluctuation integrals. The transverse 
gauge field and the Higgs field component h lead to the same contribution in both theories. 
The Goldstone channel (p fulfills a coupled differential equation. In the gauge invariant 
description we have a coupling between the field itself and its canonical momentum. In 
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the 't Hooft-Feynman-gauge, it couples to the time component of the gauge field. This 
component do not appear in the new description because they are unphysical. In the 't 
Hooft-Feynman gauge the unphysical degrees of freedom are compensated by the ghosts. 

In the same way we have found the zero mode equation, we can derive the mode functions 
for the fluctuation fields. We find 



dt 2 



u.Jt) 



u±(t) 

U h (t) 











(2.46) 
(2.47) 
(2.48) 
(2.49) 



In the mode equations for U v and , a problem in the IR region becomes obvious. The 
denominator with k 2 leads to problems for vanishing momentum. During the discussion 
of the comparison of the different approaches in section [IV], it will become clear that this 
IR-instability is caused by higher loop effects. It is possible by combining the differential 
equation for the field (|2.48|) and its conjugate momentum (|2.49|) to find an IR-stable mode 
equation by neglecting all terms of higher order than one loop. 

The energy density can easily be calculated from the Hamiltonian. We have to insert the 
expansion of the fields in dependence of the mode functions with the corresponding creators 
and annihilators. Then we have to take the expectation value in the ground state. The field 
expansion is analogous to the one for the equation of motion (|2.37 ). We find for the energy 
density in the Fourier space 



8 = -U 2 (t) + U[<p(t)] 



+ 
+ 



d 3 k e 



(2tt) 3 Aupok 2 

-n(t)<f>(t) [u: v (t)u v {t) + u^{t)u;{t)]} 



(2.50) 



It is easy to check the conservation of the energy by determining the time derivative. By 
using the equations of motion we can show that it vanishes. 
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III. COULOMB GAUGE 



A. Equation of Motion and Energy Density 

Our starting point is again the Lagrangian (|2.1p. As discussed, e.g. in M, and in 



the last section it is possible to find a gauge invariant description of the Abelian Higgs 
model by quantizing the theory with Dirac's method. An equivalent way is to choose the 
Coulomb gauge condition V • A = 0. One gets a Hamiltonian written in terms of transverse 
components and including the instantenous Coulomb interaction. This Coulomb interaction 
can be traded with a Lagrange multiplier field linearly coupled to the charge density. This 
leads to the Lagrangian in the Coulomb gauge 

£ = ^d^d^ + -8^8^ A T - eA T ■ j T - e 2 A T ■ A T &$ 

+i(VA ) 2 + e 2 ^$t$ _ ieAop _ v{&$) , 
j T = i($ty T $ - Vt^) , p = -i($$ f - $ f $) , (3.1) 

where At is the transverse component of the gauge field. The field A is a gauge invariant 
Lagrange multiplier whose equation of motion is algebraic: 

V 2 A {x,t) = p{x,t) . (3.2) 

Using the usual decomposition of $ into an expectation value and a fluctuation part, and 
splitting the fluctuations into a real part h and an imaginary part <p, we can derive the 
equations of motions for the different fluctuations and the zero mode. We give the equations 
in the Fourier space. The physical degrees of freedom - the transversal gauge mode and the 
Higgs mode h - are the same as in the gauge invariant approach and in the 't Hooft-Feynman 
gauge. We now investigate the remaining degrees of freedom, the Goldstone field p and the 
Lagrange multiplier a . As already mentioned, the equation that fixes a is purely algebraic. 
The field equation reads 

Lu 2 a (t)a (t) = e [p{t)(j){t) - <p(t)j)(tj\ . (3.3) 
For the Goldstone field we find 

<p{t) + \k 2 + A(0 2 (t) - v 2 )] ip{t) = ed o (t)0(t) + 2ea o (t)0OO . (3.4) 



It is possible to eliminate the mode function for the Lagrange multiplier a® in ( |3.4| ). By 
using the differential equation for ip and the classical equation of motion 

4> - \<j)((j) 2 - v 2 ) = , (3.5) 

it is easy to show that the time derivative of a is given by 

= W)^) - <P(tW) - 2e<f>(t)i{t)ao(tj\ = -e4>(t)<p(t) . (3.6) 

Inserting ( p73|) and (|3.6| ) in (|3.4|) leads to 

M w (t)ip(t) = , (3.7) 
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where M. w is given by 
dt 2 



+ k 2 + e z <P\t) + \ <P\t) 



2e 2 0(t) 



(3.8) 



ao is a dependent mode and the physical modes are the two transversal gauge fields a± and 
the scalar Higgs mode h. From the Lagrangian ( [3.1|) , we can also derive the equation of 
motion for the zero mode d>. We find 



<j> + \<p(<f) 2 - v 2 ) + 3\<j)(h 2 ) + 2e 2 <P(a\) + A0(<^ 2 ) 
-e 2 (p{al) + 2e(a (p) + e(d ^) : 







(3.9) 



a 



<P 



— {(f)ip - <jxp) 



where the expectation values for the fields and their normalization are not specified yet. 
Now we eliminate the field ao without using the classical equation of motion. We set <ft equal 
to 

= -A0(0 2 - v 2 ) + K(t) , (3.10) 
where TZ(t) contains the fluctuations. The field equations are then given by 

(3.11) 
(3.12) 

(3.13) 

, eX -» to r + \(<p 2 -v 2 ) . (3.14) 

By taking the fluctuation integral 7Z(t) into account, we obtain a mode function which is 
of order higher than one loop. In this case we are dealing with the back reaction of the 
quantum fluctuations ao and if. Obviously, these higher loop terms contain a factor e 2 /k 2 
as in the Hamiltonian approach which lead to the IR-instabilities. By choosing 7Z(t) = 



2e 2 



UJf, 



with 



k 2 + e 2 ^ 2 



we consider only one loop effects and ( |3.13| ) is identical with ( |3.7| ). We will discuss this case 
in detail at the end of this section. Using ( |3.11|) and ( |3.12| ), we find for the zero mode the 
equation of motion in terms of the quantum fluctuations 

<j> + X(j)((j) 2 - v 2 ) + 3X(f)(h 2 } + 2e 2 (j)(a 2 ± ) 



+ (A 



0J2 



,2ujI 



,2 A 2 



COT. 



2eV(— ^)+e^(f-) = 0. 



k 2 



(3.15) 



Again we see here the appearance of a higher loop contribution which is IR-divergent. We 
can also compute the energy density in the Coulomb gauge. It is given by 

- ,n A 



£ 



4 



(<p 2 -v 2 y 



+ 



1 r 

2 
1 



.2 2\ 



(3.16) 
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The component (a 2 ,) has a negative sign because of the indefinite metric of the time compo- 
nent of the gauge field. By inserting ( 3.11| ) into the last term of the energy, £ can be written 
similar to the equation of motion in the form 



£ 



+ 



A 



+ 



<P 



(3.17) 



With ( |3.10|) , ( |3.13|) and the equations of motion for h and a±, it is straightforward to show 
that the time derivative of £ vanishes. Now we have to decide whether we neglect lZ(t) or 
not. For a well defined one loop approximation which we have considered in the 't Hooft- 



Feynman gauge in [22] we have to take it to be zero. Numerically, this leads to some problems 
for the energy conservation. In order to show the energy conservation, we have to compute 
the time derivative of £. In the energy density terms proportional <fi appear, and, after 
performing the time derivative, it is necessary to insert the equation of motion for the zero 
mode. Since we have taken lZ(t) = 0, we have to insert <fi = — A0(0 2 — v 2 ). Analytically, this 
is no problem, but numerically this equation is only solved at the initial time. By solving 
the equation of motion for the zero mode numerically we automatically take the fluctuation 
integral into account. Therefore, we cannot expect exact energy conservation. On the other 
hand this is a good cross check for the quality of the one loop approximation. For small 
coupling constant e 2 the numerical energy conservation has to be acceptable. 



IV. EQUIVALENCE OF THE DIFFERENT APPROACHES 

At this point we can compare the gauge invariant approach and the Coulomb gauge. In 
the gauge invariant approach the Goldstone channel was described by a combination of two 
first order differential equations for the field itself ( |2.48|) and its canonical momentum ( [2.49 ). 



Now we show that this approach also leads in the one loop order to the same equations as in 
the Coulomb gauge. If we differentiate (|2.48| ) with respect to t and use the classical equation 



of motion <p = IT, we find a second order differential equation of the form 

% + (0 2 + #) U v + e ^U v - - f 2 U^ = . (4.1) 

With the relations for LL and 

U* 9 = -A + ^rU, , (4.2) 
tt e 2 00 • eVV 2 e 2 - 2 

and again with the classical equation of motion now in the form <ft = — \<p(4> 2 — v 2 ), it is 
straightforward to show 

M W U V = , (4.4) 
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where M. w is given by ( |3.8|) . By inserting the classical equation of motion without the 
fluctuation part, we suppress the higher loop terms and therefore get rid of the IR problem. 
Therefore, we have shown that if the classical equation of motion is fulfilled, the Goldstone 
mode in the Coulomb gauge and in the Hamiltonian approach has the same fluctuation 
operator. We also want to compare the equation of motion for the zero mode <j) in the 
Coulomb gauge and in the Hamiltonian approach. In order to shorten the notation and 
make a comparison easier we also write here the fluctuation integrals as expectation values 
and do not worry about the normalization. By taking the time derivative of the zero mode 
equation ( |2.39| ), we find: 



fl + rie 2 (^)+2e 2 n(^)-e 2 , 



k 2 ' 



k 2 



' k 2 



(4.5) 



Now we insert the differential equation for the conjugate momentum fl ( j2.40| ). Since we are 
only interested in the one loop order we can use the classical equation of motion = II with- 
out fluctuations and neglect terms of higher orders arising from the product of fluctuation 
integrals. Then the linearized field equation reads as a second order differential equation 



4> = -u'(<f>) - <f> 



7T 2 
2/ n ip\ 



2e 2 (ai) + 3\(h 2 ) + \(<p 2 )+e z (-t) 



e 2 f/'(0)(^ 2 ) + 2e 2 0(p-e 2 



k 2 



(4.6) 



By using ( [4.2D and (|4.3| ) for the conjugate momentum of the fluctuation field cp, we get the 
same result as in ( [3.15| ) if we choose 1Z{t) = 0. 



V. RENORMALIZATION 

We are interested in the behavior of the zero mode under the influence of the quantum 
fluctuations and in the energy. Therefore, we need well defined, finite relations, and we 
have to renormalize the equation of motion and the energy density. With a perturbative 
expansion of the mode functions it is possible to extract the leading divergences and to 
introduce counter terms independent of the initial conditions. This method was developed 
for non-equilibrium systems in [^6| and implemented for gauge theories in |p2[ . We have 
given some details of the perturbative expansion of the mode functions in appendix [A]. In 
order to make a comparison with the 't Hooft-Feynman gauge possible, we also introduce 
some renormalization conditions. We choose as renormalization point the minimum of the 
effective potential, which is due to Nielsen's theorem the same point for both gauges. The 
renormalization has to satisfy the condition |$| = v or in other words the one loop effective 
potential has to satisfy some basic properties of the classical potential. These include the 
same minimum and the same curvature at the minimum for both potentials. Therefore, we 
have to add some finite corrections to the Lagrangian. In general the effective potential is 
not very useful in the context of non-equilibrium dynamics as it is discussed for example 
T5| . Especially, the imaginary part of the effective potential leads to instabilities and to 



in 



misleading results for the investigations of the dynamics of the system. Nevertheless, the 
minimum of the effective potential point seems to be a natural point for the renormalization 
since we are using for the non-equilibrium calculations the same approximation scheme. For 
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our numerical investigations we consider only initial conditions near the minimum of the 
potential since this is the only stable point for the one loop approximation. If we would like 
to enter in the instable region of the potential then we have to take the back reaction also 
in the mode functions into account. This is beyond the scope of this paper. We are here 
interested in the comparison of different gauges and the influence of additional gauge modes 
in general. Therefore, it is useful to restrict ourselves to the stable region of the potential. 



A. Equation of Motion 

First of all, we have to write the equation of motion ( |3.15| ) in terms of the mode functions. 
Details for the quantization of the Goldstone mode ip are given in appendix [A]. In order to 
deal with a well defined one loop approximation we choose 1Z equals to zero. The equation 
of motion for the zero mode then reads 



0(t) + A cf) 2 {t)-v 2 <f>{t)+F{t) = 0, 



(5.1) 



with 
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(5.2) 



The first four integrals contain divergent contributions, which we can extract by introducing 
the perturbative expansion for the mode functions as explained in appendix Dimensional 
regularization of the divergent parts leads to the following counter terms 
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(5.6) 
(5.7) 



After introducing the counter terms we finally get the following finite equation of motion 



;i + 6Z)4> - -(m 2 h + Am 2 )0 + (A + AA)0 3 + T &n = 



(5.8) 
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with finite corrections for the wave function, the mass, and the coupling constant 
AZ = -2e 2 C , 
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The terms independent of the initial conditions are due to the choice of the renormalization 
conditions The finite fluctuation integral reads 
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Therefore, we have a well defined finite equation of motion which we can investigate numer- 
ically. 



B. Energy Density 

In order to find the divergent parts of the energy density we have to introduce the mode 
functions for the different fields in Eq. ( |3.17| ) in the same way as for the equation of motion. 
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This leads to 
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As in the equation of motion the first four integrals contain divergences. With the help of 
the truncated mode functions fh, f±, and f v we find after a lengthy but straightforward cal- 
culation in addition to the mass and coupling constant counter term a cosmological constant 
counter term which arise from a quartic divergence of the form 
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We also get an additional finite term 
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so that the complete finite energy density is given by 
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+eW (— + e \U I 2 
d 3 k 2k 2 + e 2 (b 2 



which can now be treated numerically. 



VI. NUMERICAL RESULTS 

In order to investigate the influence of the gauge field sector on the zero mode in a system 
out of equilibrium, we have carried out some numerical examples. We are interested in two 
different aspects: first we investigate the influence of the different gauges. Therefore, we 
reinvestigate the 't Hooft-Feynman background gauge which we have discussed in detail 
in P2J and we compare the results with those we find in the Coulomb gauge. Secondly, 
we investigate the effect of the different degrees of freedom which arise in gauge theories. 
For the 't Hooft-Feynman gauge we summarize the results for the equation of motion, the 
mode functions, and the energy density which we have published in ||22|| . For more details, 
especially in view of the renormalization, the reader is referred to our paper 



A. The 't Hooft-Feynman Gauge 

We give a short review on the main steps to derive the equation of motion and the 



energy in the 't Hooft-Feynman background gauge. In opposite to our previous work [22 
we restrict ourselves here to the abelian Higgs model. This leads to some modifications in 
the degeneracy factors which we have introduced in ]22|. The Lagrangian which describes 



the abelian Higgs model in the 't Hooft Feynman gauge contains three parts 

Aot = £ + £gf + £ F p ■ (6.1) 
£ is given again by ( j2.1|) . We split the Higgs field into a mean value and fluctuations 

$(£, t) = [<f>(t) + h(x, t) + i(p(x, t)} , (6.2) 
and consider for the gauge field only fluctuations 

An(x,t) = a M (f,t) . (6.3) 
The gauge fixing term has the following form 

C GF = -\f 2 , (6.4) 

with 

% lV )=a/ + e(^ + %. (6.5) 
The corresponding Faddeev-Popov Lagrangian which is relevant for our calculations is 

£ FP = d^rfdr] — e 2 4> 2 rfrj . (6.6) 
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After inserting the expansions for the fields in the total Lagrangian (6.1) and neglecting 
all terms higher than second order in the fluctuations we derive the following equation of 
motion for the zero mode 



+ A0(0 2 - v 2 ) + 3\<f)(h 2 ) + e 2 <f)(a 2 ± ) 

+ (A + e 2 ) (/>(<p 2 ) - e 2 <J)(a 2 ) - ed t (a ip) = . 



(6.7) 



In this notation we have not taken care of the normalization and the renormalization or 
the different solutions for the mode functions of the coupled channel. We only want to give 
an overview of the equations and the connections of the fields and do not go into technical 
details. As already mentioned, the mode functions for the isoscalar Higgs field and the 
transversal gauge field are the same as in the Coulomb gauge ( |2.46| ) , ( |2.47| ) . For the coupled 
sector ao(p, they are given by 



g<P(t) 




. 



(6.8) 



The operator has two interesting features which distinguishes it from the single modes, the 
indefinite metric of the time component of the gauge field, and the time derivative of the 
zero mode in the off diagonal elements which connect the two fields. We find analogous 
properties in the Coulomb gauge; the fluctuation part for ao contributes with a negative 
sign to the fluctuation integral ( |3.9| ), and we find time derivatives of the zero mode in the 
equation of motion for the zero mode ( [3.15D as well as in the mode function for (p ( p.7| ). 

In the Feynman gauge, the mode functions for the transversal gauge field, for the longitu- 
dinal gauge field, and for the ghost fields are the same. Two of the three gauge components 
are canceled by the ghost fields and only one degree of freedom is left in contrast to the 
Coulomb gauge where we have two transverse gauge components. The energy density can be 
derived by integration of the equation of motion for the zero mode or from the corresponding 
Hamiltonian of the system. It reads (see also 
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1 12 , X 
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(6.9) 
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1 r„- 



+ - 



(0 2 > + KV> 



K> + {<<)_ 

(«0> + (Vl<4) 



This expression looks very similar to the Coulomb energy fl3.16| ) despite the fact that ao is 
dynamical and contributes with a derivate part. 



B. Results 

For our numerical calculations, we have chosen four different sets of parameters listed 
in Table Q. The initial value for the zero mode <ft an d the Higgs mass are the same for 
all sets. They are chosen in such a way that the zero mode evolves in the right minimum 
of the potential. With this choice of initial conditions, the zero mode can not evolve into 
the complex part of the effective potential around zero. In this region, instabilities increase 
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dramatically and the one loop approximation breaks down as explained for the 4 theory 



for example in [ 15 . [2! 

In order to give an impression of the potential we have plotted in Fig. |I] for the Coulomb 
gauge and Fig. for the Feynman gauge different approximations for the potential. The 
dashed line shows the classical potential V((p) = ^(</> 2 — v 2 ) 2 , the circles include finite 
corrections due to the renormalization, and the pluses display the effective one loop potential. 
Since we have chosen the renormalization conditions in such a way, that the classical and 
the effective potential have the same minimum and the same curvature at the minimum, 
they are very close to each other. The solid line shows the zero mode part of the energy 
versus <p. The field begins to roll down the potential but the energy is not high enough for 
the field to reach the maximum at zero. Therefore, it starts to oscillate in the minimum. 
The prediction of these two plots in the context of non-equilibrium dynamics is not clear, 
because the effective potential is an equilibrium quantity, since the expectation value of the 
scalar field, that serves as order parameter, is space time independent. Nevertheless, they 
are instructive to get an idea of the potential by which the zero mode is influenced. 

For our numerical considerations we have only varied the coupling constants A and e and 
therefore the masses of the different fields. We have also given the initial masses of the three 
different fields in Table |. 

For the first parameter set, we have chosen the same coupling constant for the Higgs field 
and the gauge field. The initial masses for the fields are all small but not zero. Since we have 
taken the initial value for the zero mode to be small, the effect of the quantum fluctuations 
is negligible and therefore, oscillates with constant amplitude. The behavior of the zero 
mode, which is displayed in Fig. [3] agrees in both gauges excellent. 

In the second parameter set we have chosen a dominant gauge coupling, therefore the 
initial mass of the Higgs field becomes very small compared to m 2 and ml xo . In this case, 
the zero mode is moderately damped, as shown in Fig. the behavior is comparable to a 
purely scalar theory. The two gauges agree again very well, at late time the zero modes are 
slightly dephased. 

The situation changes drastically, if we choose a smaller gauge coupling and therefore 
a nearly vanishing Goldstone mass. In Fig. |5| and Fig. |6] we show two examples for this 
constellation. 

In the first example we have chosen a nearly vanishing gauge coupling with e = 0.001. 
In this case an instability occurs and the system breaks down. We have analyzed a similar 
problem occurring in the 4 theory in |HJ. There we have found that the spontaneous 
symmetry breaking leads to instabilities if the mass squared term of the fluctuation field 
becomes negative due to the decrease of the zero mode. In addition to m1(t) = ml + 
3A(0 2 — v 2 ), the mass term, which is also relevant for the spontaneously broken 4 theory, 
we now have to secure, that rn 2 x (t) = e 2 <p 2 + A(0 2 — v 2 ) does not become negative when (p(t) 
decreases. This leads to the two conditions 

m 2 h (t)>0 =► <P(t)>^, (6.10) 

2 

m 2 eX (t) > =► <j>(t) > . 2(e2 * A) . (6.11) 

Therefore, the critical value for the field <f>(t) for parameter set 3 is around (p(t) > 0.5 (exactly 
0.49999975) in order to find a stable development for the system. We have displayed in the 
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frame the development of the zero mode for a shorter time period. The solid line at 0.5 
shows the critical value for 0, for which rn 2 eX (t) becomes negative. One sees that the average 
value of the zero mode is slightly smaller than this critical value. Therefore, the system 
breaks down. 

In the second example we have chosen the gauge coupling to be e = 0.1, therefore slightly 
larger but again small. In this case the zero mode is strongly damped and settles down to 
the minimum. An appropriate size for the gauge coupling prevents the system to destabilize. 
We can influence m e \(t) either by changing the Higgs coupling or the gauge coupling. As 
this example in Fig. ^| shows it is possible to find a proper relation between A and e in order 
to allow the decay of h into (p and gauge field and a stable development of the system in the 
minimum of the potential. We have also plotted the critical value for the zero mode in this 
case where it is obviously much smaller than the final value of the zero mode. 

In order to understand this strong damping of the field more precisely, it is instructive to 
investigate the Feynman graphs for the problem in detail. This is especially easy in the 't 
Hooft-Feynman gauge. We have listed the the occurring graphs in appendix [TJ. The relevant 
graph which leads to the strong damping is the following 

<p 

h o 



El 

It allows the decay of the Higgs field into a and <p for a proper relation between the 
masses. In the Coulomb gauge this mechanism is not as obvious as in the 't Hooft-Feynman 
gauge since we have eliminated the ao-mode. Also the structure of the fluctuation integral 
for the p-mode is rather complicated and allows therefore not such a simple analysis. 

In order to check our numerics we have plotted the energy density for the Coulomb gauge 
in Fig. |7| and the 't Hooft Feynman gauge in Fig. § for the fourth parameter set. The upper 
line displays the fluctuation energy which increases and the lower line the zero mode part of 
the energy which decreases. The solid line shows the total energy. For convenience we have 
added in both cases a constant to the fluctuation part of the energy. Otherwise the curves 
are only straight lines due to their distance. The energy conservation is excellent in both 
cases. 

Summarizing the results, we have found that the damping effect is strongest for a nearly 
massless Goldstone field and a massive isoscalar Higgs field. In this case the isoscalar Higgs 
field has the possibility to decay into the other fields. We have found an analogous behavior 
in the 4 theory in the large N limit. There, the damping of zero mode was due to the 
massless Goldstone bosons. 

We have also found that the behavior of the zero mode and the energy is the same for 
the different gauges and that the Goldstone channel plays an important role. 

VII. CONCLUSIONS AND OUTLOOK 

We have analyzed in this paper the abelian Higgs model out of equilibrium in detail. 
Therefore, we have investigated a gauge invariant approach which was developed by Boy- 
anovsky et al. [[TJ]. We have extended their calculations to a complete set of equations, which 
describes the evolution of the zero mode under the influence of gauge and Higgs fluctuations. 
We found some problems in the IR region induced by the inclusion of terms higher then one 
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loop order. Since these terms were included in an uncontrollable manner, we have computed 
a linearized form of the equations. They are equivalent to the Coulomb gauge fixed theory 
in the one loop approximation. We have performed the renormalization of the system in 
order to get finite well defined equations which we have investigated numerically. We were 
also interested in the behavior of the system under the influence of different gauges. Based 
on [22 1 we have reinvestigated the 't Hooft-Feynman gauge. The main differences between 
the two gauges is that we have eliminated in the Coulomb gauge all unphysical degrees of 
freedom while in the 't Hooft-Feynman gauge we add extra degrees of freedoms, the ghosts 
fields, in order to cancel the unphysical degrees of freedom. We found a very good agreement 
in both gauges for the development of the zero mode and the energy. The influence of the 
various degrees of freedoms which arise in gauge theories have led to new and interesting 
results in the behavior of the zero mode. Due to the possibility for the Higgs field to decay 
into other fields, the zero mode was efficiently damped. This behavior was until now not 
observed for systems treated in the one loop approximation beside for the case of fermion 
decay p7 |. As in the one loop approximation in the </> 4 theory the evolution of the zero 
mode in the unstable region was impossible. 

For future studies there are two different interesting aspects. The first one is more tech- 
nical. In order to enter the unstable regime of the effective potential it is necessary to 
go beyond the one loop approximation. Therefore, the implementation of Hartree-like ap- 
proximation would be very interesting. As we have shown here one has to include higher 
order effects in a careful way in order to avoid IR-problems. In addition one has to choose 
very carefully an expansion parameter in order to secure gauge invariance. Since the usual 
Hartree-approximation is just an ansatz and an inconsistent summation and not an ex- 
pansion in the coupling this will probably lead to problems. Also the investigation of the 
non-abelian gauge would be very interesting. Beside these technical considerations, the im- 
plementation of our results for the gauge fields in a cosmological context would be very 
interesting. Together with our studies on fermionic systems out of equilibrium we now have 
built the foundation to examine more realistic models to describe the physics of the early 
universe. The recent success in detecting neutrino masses has revived the idea of grand 
unification. An implementation of our method in Grand Unified Theories could lead to new 
and interesting results. Also the implementation of Friedman- Robert son- Walker cosmology 
in the model we have considered is interesting in order to get a more suitable model for 
describing the inflationary scenario. 
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APPENDIX A: PERTURB ATIVE EXPANSION 



In this appendix we give some details of the perturbative expansion of the mode functions 
which we have used for the renormalization procedure. We have developed this formalism 
for non-equilibrium systems in [^BJ and extended it for gauge theories in P2 |. For more 
details the reader is referred to these works. First of all we have to quantize the Goldstone 
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field (p and introduce the mode representation in the same way as the transversal gauge 
mode and the Higgs mode in Eqs. (|2.46|) and ( 2.47 ). Therefore, we have to investigate the 
following equation for the Goldstone field fl3.7| ) 

<p + ul x (t)<p-^\<l><p-j><p\ =0. (Al) 

In order to quantize the field we have to find the conjugate momentum. We can read it of 
from the Lagrangian fl3.1|) : 

, k 2 e 2 (j)(j) . . . 

11^ = y? - ea (f) = tp— —Lp . (A2) 

The commutation relation for the field and its momentum is given by 

[<p,H v ] =i5(x-x > ) . (A3) 

By computation of the commutator for the field and its time derivative, we find that it is 
multiplied by a factor 

k 2 

[<p,<p]— = i5(x-x'). (A4) 

Now we can expand the field in terms of the mode functions and the corresponding annihi- 
lation and creation operators 

/ii 3 k 1 - - 

(2^8 ^ ("Mt)e a * + a\u;{t)e-^) . (A5) 

In order to satisfy the commutator for p> and tp, the commutator for a k and a\ differs from 
its common form 

[a k ,al] = 2i(2n) 3 5 3 (k-k')^. (A6) 
The expectation value for the field is then given by 

w = /Sw- (A7) 

The mode function satisfies the differential equation 

U 9 + lu 2 x U v - ^ [<j>U v - j>U v ) = . (A8) 

The frequency in the expectation value has to be fixed by the determination of the Wronskian 
which belongs to the mode equation. Since the form of the differential equation is not 
standard due to the first derivative of U v , the following transformation is efficient 

U v = uj a U v . (A9) 

The f/^-terms are canceled and we find a new mode equation of the form 

~ / „ 3e 2 6 2 k 2 e 2 M\ ~ , . 

U v + u? eX + — ^— + )U V = . ( A 1 ) 
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Now we redefine the frequency uj 2 x in order to find a suitable Wronskian. We choose 



_ 2 _ 2 3eW e 2 # 
and get the new differential equation for the mode function 

The Wronskian then has the well known form 

6 v u; - uJj\ = c . 



(All) 



(A12) 



(A13) 



Since U v behaves like e ^ eA0 * we can fix the constant to be C = —2iu e \ . For the initial 
time, and the fluctuation integral vanish, so that u 2 X0 simplifies to 



u 2 



e 2 X 



exo = k + z% + Kn - v) - —n 



(A14) 



For high momenta, which are important for the UV-divergences, the last term is negligible. 
Now all mode equations have the similar structure 



d 2 2 

dp + Wj0 



dt 2 



+ ^eAO 



Uj(t) = -V 3 {t)U 3 {t) , 

U v (t) = -Vex{t)U^{t) , 



with j = a, h. Here we have introduced the potentials 

V a (t) = e 2 [0 2 (t)-02] ; 
V h (t) = 3A [0 2 (t) - <j>l] , 

yr (f] v (f] , Se 2 ^t) 2 k 2 , e 2 0(t)0(t) e 2 A0 2 (0 2 - v 2 ) 

V e \\t) = V e \{t) H 7777 1 7777 h 



V;a = (e 2 + A) (0 2 - 2 ) 
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(A15) 
(A16) 

(A17) 
(A18) 

(A19) 
(A20) 



We separate the mode functions into a trivial part corresponding to the case that the poten- 
tial vanishes and a function fj(t) which represents the reaction to the potential by making 
the ansatz 



Uj(t) = e-*>»*[l + fj(t)} , 
%{t) = e'^ [1 + U(t)} . 

Then the new mode functions satisfy the differential equation 

f j (t)-2kj S oMt) = -V S (t)[l + fM, 



(A21) 
(A22) 

(A23) 



respectively for with the initial conditions fj(0) = fj(0) = 0. Expanding now fj(t) with 
respect to orders in the potential by writing 



m=fr+fP+fP + 
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(A24) 



we can extract the leading behavior for fj(t): 

ffit) = ~Jdlf W) -W + 4^J dt> ^^V) + 0{^) , (A25) 
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with At = t — t'. For f v the frequency and the potential have to be replaced by u e \ and 

VeX. 

Inserting the leading behavior of the new mode functions into the fluctuation integral 
( [5. 13] ) allows us to extract the following divergences 
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Dimensional regularization leads to the counter terms given in Eqs. ( p.3| )-( p75|) and finite 
corrections. 

In a similar way we have to treat the energy density ( |5.14 ). Here we find with the help 
of the truncated mode functions / 
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w eA0 , VeX | ™a™?> 



4cj 



eAO 



4A; 2 



^eAO 



16^ e 3 A0 



4u^ n £> 



aO^eAO 



which leads again to Sm, SX, and SZ and in addition to 5A given in Eq. (|5.15 ) and also to 
an extra finite contribution. 



APPENDIX B: THE LEADING FEYNMAN DIAGRAMS 

In order to understand the behavior of the system it is very instructive to analyse the 
leading Feynman diagrams in detail. This is especially easy in the 't Hooft-Feynman gauge, 
since we can read of the Feynman rules directly from the gauge fixed Lagrangian, which is 
given for example in | 2~2]| . We find four different propagators which are 
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1. the gauge boson propagator: 



iA 



ab _ ig^S ab 



"w 

with m 2 ^ = e 2 v 2 , 



(Bl) 



2. the propagator for the isoscalar Higgs field: 



with ml = 2Xv 2 . 



(B2) 



3. the propagator for the Goldstone field: 



iA ab - i K 



(B3) 



4. the propagator for the ghost field: 



(B4) 



and ten vertices which are relevant if we restrict ourselves to second order in the fluctuations 
VI) V2) 





*r = -z|A [0 2 (t) - v 2 



V3) 



-i (| + f ) 5 ab l<p 2 {t) - v 2 



V4) 




7f 



iT = -ie 2 5 ab [0 2 (t) -v 2 ] 
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V5) 



V6) 




iT = ie 2 e abc q^{t) iF = 2ie5 ab q fl <f){t) 

These lead to the following Feynman graphs up to second order in perturbation theory 

± +[ (B5) 
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- - 

-V- 



The plus and minus signs at the graphs indicates the use of propagators in the CTP- 
formalism. Especially interesting is the last graph; it leads to the possibility of the decay 
of h into ao and (p which causes the strong damping of the zero mode for parameter set 
4 displayed in Fig. ||. The various degrees of freedom in gauge fields lead to new and 
interesting results in non-equilibrium quantum field theory. A detailed analysis of these 
graphs within the CTP-formalism is given in p4 . 
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A 


e 2 


m\ 


0o 






m lxo 


Parameter set 1 


1 


1 


0.5 


0.51 


0.26 


0.53 


0.27 


Parameter set 2 


0.33 


1.3 


0.5 


0.51 


0.44 


0.01 


0.28 


Parameter set 3 


1 


0.001 


0.5 


0.51 


2.6-10" 7 


0.53 


0.01 


Parameter set 4 


1 


0.1 


0.5 


0.51 


2.6-10" 3 


0.53 


0.01 



TABLE I: Parameter sets 
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FIG. 1: Potential vs. <p in the Coulomb gauge for parameter set 4, solid line: classical energy vs. (f>, 
dashed line: classical potential, circles: classical potential with finite corrections, pluses: effective 
potential. 
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FIG. 2: Potential vs. (f> in the Feynman gauge for parameter set 4, solid line: classical energy vs. <fi, 
dashed line: classical potential, circles: classical potential with finite corrections, pluses: effective 
potential. 
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FIG. 3: Zero mode vs. t for parameter set 1, dashed line: Feynman gauge, dotted line: Coulomb 
gauge. 
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FIG. 4: Zero mode vs. t for parameter set 2, dashed line: Feynman gauge, dotted line: Coulomb 
gauge. 



30 




FIG. 5: Zero mode vs. t for parameter set 3, dashed line: Feynman gauge, dotted line: Coulomb 
gauge; frame: zoom onto a shorter time scale, solid line: critical value for 4>{t). 
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FIG. 6: Zero mode vs. t for parameter set 4, dashed line: Feynman gauge, dotted line: Coulomb 
gauge, solid line: critical value for (f)(t). 
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FIG. 7: Energy vs. i in the Coulomb gauge for parameter set 4, solid line: total energy, dashed 
line: classical energy, dotted line: fluctuation energy. 
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FIG. 8: Energy vs. t in the Feynman gauge for parameter set 4, solid line: total energy, dashed 
line: classical energy, dotted line: fluctuation energy. 



32 



